Out of equilibrium thermal Casimir effect in a model polarizable material 
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i-H ' 

■») ' Relaxation of the thermal Casimir or van der Waals force for a model dielectric medium is inves- 

, tigated. We start with a model of interacting polarization fields with a dynamics that leads to a 

frequency dependent dielectric constant of the Debye form. In the static limit the usual zero fre- 
quency Matsubara mode component of the Casimir force is recovered. We then consider the out of 
I ^ | equilibrium relaxation of the van der Waals force to its equilibrium value when two initially uncor- 

related dielectric bodies are brought into sudden proximity. It is found that the spatial dependence 
of the out of equilibrium force is the same as the equilibrium one but it has a time dependent ampli- 
tude, or Hamaker coefficient, which increases in time to its equilibrium value. The final relaxation 
f-H , to the equilibrium value is exponential in systems with a single or finite number of polarization field 

relaxation times. However, in systems, such as those described by the Havriliak-Negami dielectric 
' constant, with a broad distribution of relaxation times, we observe a much slower power law decay 

to the equilibrium value. 

+2 . I. INTRODUCTION 

Despite the success of theoretical calculations on the equilibrium Casimir force [H-Q , aspects of the out of equilibrium 
behavior are still poorly understood and the subject of controversy and debate. A number of approaches have been 
adopted to compute thermal fluctuation induced forces out of equilibrium in simple models of soft matter systems and 
binary liquids . For example, the stress tensor has been used to compute the force [34H in a variety of non-equilibrium 
contexts. While it is clear that in such systems computations using the stress tensor will yield the average value of 
the force at thermal equilibrium, it seems nevertheless that more information is needed regarding the dynamics of the 
field theory representing the critical or fluctuating field |9( , in particular how the value of the field at a surface changes 
when the surface is moved. An alternative approach is to define an energetic interaction of the field with a surface 
and then define forces via the principle of virtual work [lfj El . Yet another is to define the force at a surface by a 
local kinetic arguments, for instance by using the ideal gas form for the pressure as a function of the local density field 
[l2l [l3j . There are notable differences in out of equilibrium forces computed using the approaches above and relatively 
, few systems have been studied explicitly. However studies of free Gaussian field theories undergoing model A (non 
conserved) dynamics have been carried out. Model A dynamics for the field is basically a diffusion equation driven 
by white noise and the thermal Casimir forces are found to tend toward the equilibrium value with diffusive scaling 
[fl S E3, El • As well as studying the approach to equilibrium for dynamics obeying detailed balance, one can examine 
£SJ | what happens when the noise is non-thermal, for instance colored. In this case the steady state Casimir interaction, 
under model A type dynamics, tends to acquire an additional screening due to temporal correlations in the noise 
0) E3) El • We also note that the Parisi- Wu stochastic quantization scheme can be used to extract equilibrium results 
for the quantum Casimir force via a Langevin dynamics approach |14j . 

In the case of quantum electrodynamics the definition of the instantaneous force can be derived from the stress tensor 
as it is physically defined, within the theory of electromagnetism, via the force acting on charges and currents just as 
was done in the first Lifshitz paper on the equilibrium van der Waals force [15| . The study of the out of equilibrium 
quantum electromagnetic Casimir effect is however much more complicated than the equilibrium thermal Casimir 
effect and a number of physical assumptions need to be made if one wants to study dynamics or non-equilibrium 
situations. One approach is to use generic models of fluctuating electrodynamics (or stochastic Maxwell equations) a 
la Rytov [l|| , where the field is driven by randomly fluctuating current densities or dipole fields [l7i - fT9l l2ll . [22| . The 
out of equilibrium context analyzed in the above studies is where the interacting media have different temperatures. 
In some cases there can be a long range bulk induced interaction between the media, essentially due to the difference 
in their blackbody radiation. A similar bulk effect was found in 5] for the thermal Casimir effect in fluctuating 
scalar fields in the presence of temperature gradients. While the theory of fluctuating electrodynamics is very general 
and applies to materials of generic dielectric properties, there are nevertheless certain local equilibrium/fluctuation 
dissipation properties that need to be assumed in these theories. The goal of this paper is not to study non-equilibrium 
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steady states but rather to address the question of how the equilibrium thermal Casimir force evolves in time to its 
equilibrium values, for instance when two dielectric slabs are brought suddenly into proximity. 

Specifically we will examine the out of equilibrium behavior of the thermal part of the electrodynamic Casimir 
force (corresponding to thermal van der Waals interactions). Our approach will be based on a specialized form of 
fluctuating electrodynamics based on the full n-body dipolar interaction between two model dielectric media. The 
question we will address in this paper is how the thermal van der Waals interaction between two media relaxes to 
its equilibrium value? Some physical properties of the thermal Casimir interaction look less mysterious when studied 
in this setting and we will show that the thermal Casimir interaction is induced by the correlations between the 
polarization fields of the two media. For a single type of polarization field our model corresponds to a Debye type 
dielectric material. However, if we assume that there is a distribution of polarization fields with different relaxation 
times and dielectric susceptibilities, then any dielectric function can be obtained by taking a suitable distribution. 
Thus although our model only applies to Debye-like dielectrics it can be applied to any dielectric function as results 
of our computations only depend of the frequency dependent dielectric function. Furthermore as we are interested in 
the approach to the equilibrium force we argue that the slowest relaxing polarizability fields will be the thermal ones 
(zero Matsubara frequency term) and among these thermal modes those with the slowest dynamics should exhibit a 
Debye relaxation. 

The dynamics we assume for the microscopic model for a polarizable dielectric media is over-damped stochastic 
dynamics for the polarization field. In the static limit, the force between such media corresponds to the zero frequency 
Matsubara term in the interaction energy between dielectric media as found in the Lifshitz theory [1[ . We show here 
how the thermal van der Waals force between two objects emerges via the correlation of dipoles in the interacting 
media. Within this formalism we can see how the force evolves with time towards its equilibrium value, for instance, 
for two slabs brought into close proximity where initially they had infinite separation. We can also see how on upon 
changing the temperature of the system, the van der Waals force evolves from its initial equilibrium value at the initial 
temperature to the final equilibrium value at a different temperature. Although the model is simple, and restricted to 
the non-quantum part of the van der Waals interaction it may give useful intuition about fluctuation induced forces 
out of equilibrium. An appealing aspect of this approach is that the Laplace transform (with respect to time) of 
the dynamical force can be expressed in terms of static results that can be obtained via well established equilibrium 
methods. 



II. MODEL OF POLARIZABLE MEDIA 



In this section we define a simple model for the polarizable media and their dielectric properties, showing that 
the dielectric response functions are a linear combination of Debye-like terms. We then show how the average force 
between two such media in thermal equilibrium is identical to the thermal Casimir or van der Waals force as predicted 
by the Lifshitz theory. 

Consider a system with an ensemble of local polarization fields Pi/(x) at the point x in the medium. The index v 
corresponds to a type or species of polarization field which notably has its own polarizability per unit volume denoted 
by x^( x )- We consider two semi-infinite regions (slabs) V + and V~ defined via the sign of the coordinate z, such 
that z > in V + and z < in V~ . The two regions V + and V~ will be separated in the z direction by a distance L. 
In this notation the total energy for a given configuration of the dipole fields is 

H = \ I d x d y p " W ' A uu , (x, y, L) Pl/ , (y) (1) 

vv' 

where 

A„„(x,y) = ^ (x ~ y)J ^' +J>(x,y,L), (2) 

In the first term / is the identity matrix in ZD space and the polarization energy Eq. Q] corresponds to the classical 
harmonic energy needed to generate a local polarization field. In the second term D is the interaction between the 
dipoles in two semi-infinite regions which we will write in terms of the standard dipole-dipole coupling 

D (x-y) = -V i V i G(x-y), (3) 

where G is the Green's function for the electrostatic field in vacuo obeying 



e V 2 G(x) = -5(x). 



(4) 
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Writing the energy in the above manner means that the separation (L) dependent part of the energy is encoded in 
the dipolar interaction term D. This interaction is given by 

£>(x,y,L)=D (x-y), (5) 

when x G V + and y G V + or when x G V~ and y G V" - , 

I>(x,y,L) =J D (x-y-Lz) (6) 

when x G V - and y G and 

I>(x,y,L) =J Do(x + Lz-y) (7) 

when x G V + and y G V~ . 

Before proceeding with the analysis of Casimir interactions between two slabs, let us consider the dielectric properties 
of the bulk systems. We will assume that each dipole field obeys Langevin dynamics and that the only coupling between 
the different polarization fields are their mutual dipolc-dipole interactions. The dynamical equations for each field in 
the presence of a uniform time dependent electric field E is thus given by 

^^ = -n u ^)-^L-+ Xu E l {t)+U(^t) (8) 
at op iv (pc) 

where k v is a local diffusion constant for each polarization field and determines the time-scale of relaxation of the 
field. The noise (x, t) is of a white noise type and has a space-time correlation function 

(Crf(x,t)£^-(x',t / )) - 2T^(x)V(x,x'), (9) 

where T is the temperature (imposed on the polarization field by the local bulk environment) and the weighting 
Kv assures that detailed balance holds, i.e. that the dynamics will eventually lead to thermal equilibrium at the 
temperature T. 

In a bulk system (where k v and \v are constant) we determine the dielectric properties of the model by computing 
the response to a spatially uniform applied electric field E(t) = Ecxp(iwt). The average response of each polarizability 
field is uniform in space and the interactions between dipoles averages to zero. The average value of each polarizability 
field is then given by 

(p„(t)) = Xh u — EexpM), (10) 

and the average total electric displacement is thus 

D(t)=e E(t) + £>„(i)>. (11) 

This allows us to read off the frequency dependent dielectric response as 

£ H = e « + E T ^7' ( 12 ) 

which is obviously just a superposition of Debye-like dielectric responses. The frequency dependent dielectric response 
thus depends on both the polarizabilities \v and the dynamical variable n v . The static dielectric constant however 
depends only on the static polarizabilities 

e(0) = e + ]Tx„. (13) 



III. EQUILIBRIUM VAN DER WAALS INTERACTION FOR SLAB GEOMETRIES 



After having established the basic bulk dielectric properties of our model we will show how they enter the equilibrium 
van der Waals or thermal Casimir effect between two semi-infinite dielectric slabs. 
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The force between two semi-infinite slabs separated by a distance I for any configuration of the polarizability fields 
in the two slabs is given by 



(14) 



since only the dipole interaction term D depends on L. The local polarizability only depends on the coordinates 
within the two volumes V + and V~ . The equilibrium value for the average force may be obtained using the fact that 

(p„(x)jv(y)) =T[A-\^y,L)U, (15) 

and thus 

(/> = -^Tr (a~^a), (16) 



2 V dL 

where Tr indicates the trace over the operator, the spatial and the species indices. This average force may thus be 
written as 

t f) fit? 
lf}= -.-T lblA] = .- (17) 

which of course agrees with the standard statistical mechanical definition if F = — Tln(Z) is interpreted as the free 
energy obtained from the partition function 

Z = J <2[p]exp(-/3if). (18) 

This partition function can be written in a standard way by introducing a Hubbard-Stratonovich auxiliary field </>, 
which physically can be identified with iip where ip is the fluctuating electrostatic potential, to decouple the dipolar 
interaction. This gives up to a constant factor independent of L 



Z = 



P,(x£ 
(x) 



J] d[p„M0 exp ( -t [ e [V<Xx)] 2 + P [ dx^ p(x) • V0(x) - V ^ 



We note that the variables \v are n °t necessarily the same in the regions z < and z > L as in general they 
will correspond to two different materials. As the integrals over p„ are now decoupled, they may be carried out to 
yield, again up to constants independent of L, the standard form of the partition function for the thermal Casimir 
contribution for van der Waals interactions between dielectric media 



Z = J d[</>] exp Jdxe(x, 



O)[V0(x)] 2 . (20) 



Here e(x, 0) = eo + x( x ) f° r z < and z > L i.e. in the volumes V + and V" and e(x) = eo for z £ [0, L] i.e. the 
vacuum between the two media. We also see that the variables e(x, 0) are simply the local static dielectric constants 
as defined by Eq. (|13[) . This form of the partition function corresponds to the one stemming from the Lifshitz theory 
of the Casimir force for the thermal van der Waals component of the interaction. 

The partition function in Eq. (|20|) can be evaluated exactly for systems where e(x) depends only on the coordinate 
z (slab-like configurations), while for other geometries it can be evaluated via different approximation schemes such 
as the proximity force approximation or using systematic multipolar expansions. 

IV. DYNAMICS 

We now turn to the problem of the dynamical evolution of the force. In the absence of an applied field the 
polarization dynamics in the two slabs can be written as 



(21) 
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where k„(x) determines the local relaxation in the region x of space. The condition of detailed balance implies that 
the noise correlator obeys 

(C iw (x ) t)C J v(x',0) =2TJi i „ J v(x,x') (22) 

where the operator R is given by 

Rivjv'fa, x') = K„(x)<5y(5j,„'<5(x - x'). (23) 
In operator notation the dynamical equations can be written as 

3p(x) 



dt 



-JL4p(x) + C(x,t). (24) 

The average value of the dynamical force Eq. (TT4|) can be obtained from the time correlation function of the dipole 
field defined as 

(pj„(x,t)p,,v(x',*)) = C w jv'(x,x',£) (25) 
so that the time dependent average force is given by 

(f(t,L,T)) = ~TrC(t)-^A. (26) 

The evolution equation for p, being of first order in time, can be integrated to give an explicit form for the correlation 
function 

C(t) = exp(-tRA)C(0) exp(-tAR) + TA _1 (1 - exp(-2tAR)) (27) 

in operator notation, where C(0) is the value of the correlation function at t = 0. We notice here that the force 
between the two regions depends on the cross correlation between them. If at t = the two regions are brought into 
proximity from a large distance the cross correlation at t = is zero and only the second term in Eq. (1271) remains. 
In the absence of initial correlations the Laplace transform of C 



CC{s) = / dt C(t) exp(-si) (28) 
Jo 

is given by 

CC(s) = -[A+^]-\ (29) 
s 2 

This means that the Laplace transform of the time dependent average force can be written as 

£</>(.) = -^+f^l]-^A (30) 

Two things should be noted at this stage: (i) the operator R does not depend on the distance between the two regions 
V + and V~ and (ii) its inverse is simply 

irWW = ^r 5 ( x - x ')- ( 31 ) 

The first of these points means that we can write 

C(f)( S ) = --^ = -^-HZ s ) (32) 
s oL s oL 



where 



Zs= jY[d[p,}ex P (~/3H d ( s )) (33) 
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and Hd(s) is an effective dynamical Hamiltonian given by 

4 ^ K„(x) ' 

It can be written as a static Hamiltonian with local dynamical polarizabilities 

1 1 s 



X«(x,s) 2k„(x) 
which in turn leads to local dynamic dielectric constants given by 



(35) 



ed (x, s) = eo+ £-^L_. (36) 

The observant reader will immediately recognize a similarity between the form of Eq. (|36l) for the dynamical dielectric 
constant and the frequency dependent dielectric constant predicted by the dielectric response model from Eq. (|12p . 
indeed we find that erf(s) = e (—if)- Thus for the computation of the thermal van der Waals forces in this model, 
knowledge of the frequency dependent dielectric constants allows one to predict the temporal evolution of the force 
towards its equilibrium value. Although the computation above was carried out with slab geometries in mind it is 
easy to see that it applies for general geometries where dielectric objects are separated by vacuum. This means that 
we can write the time dependent force as 

(/(*)> (37) 

where F (t, L) is an effective time dependent free energy given by 

F{t,L) = -TL- 1 - s \n(z{e{-i S -))), (38) 

where L^ 1 indicates the inverse Laplace transform and the notation e(is/2) denotes that the frequency dependent 
dielectric constant is taken in all regions at the value —is/2. The pole at s = yields the equilibrium free energy and 
thus the equilibrium force. In addition we can invert the Laplace transform using the Bromwich integration formula 
to get 

F(t, L) = -T £° J^- exp(st) In (z(e(-t|)) . (39) 

where the integration is to the right of the imaginary axis (as all singularities arc at negative s). One can remove the 
term with the pole at s — by hand to write 

F(t, L) = F eq (L) - T £2^- s exp(si) [in - ln(Z(e(0))] (40) 

where 

F eq (L) = -T\n(Z(e(0)) (41) 

is the equilibrium free energy F eq (L) = lim t _ ) . C)0 F(t, L). In the above formula the remaining contour integral is 
now free of singularities on the imaginary axis and we may therefore write it as a Fourier transform by making the 
substitution s = 2iu 

-^exp(2iwt)K^(eH) - HZ(e(0))}. (42) 

- oo 27rt w 

The above results can be generalized to the situation where the two slabs are in equilibrium (at fixed distance L) 
at a temperature To and the temperature is then changed to T. We can write the force for two slabs as the sum of 
two components 

(f(L,t,T -> T)) = (f(L,T )) eq + {T T T °\ f^(t,L,T)) (43) 

where (f(L,To)) eq is the equilibrium force at the temperature To and (f^(L,t,T)) is the time dependent force for 
two initially uncorrelated slabs at temperature T. The result for uncorrelated slabs at temperature T can obviously 
be extracted from the above result by setting To = 0. 
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V. ANALYTICAL RESULTS FOR SLAB GEOMETRIES 



Here we consider the case of two parallel semi-infinite slabs of dielectric constants e\ and ti (we take the subscripts 
1 and 2 to refer to the regions V~ and V + respectively). The dynamical free energy 

F s = -T\n(z{e{-i S -)) (44) 

can then be read off from standard equilibrium results and is given by 



TS f°° / s s 

Fs= 16^ J Udu ln ( 1_ Al (- J 2 )A2(_J 2 )eXp( ~ u) 



(45) 



where S is the area of the slabs and 



A l(w ) = e f\- e ° (46) 

The time dependent force for two initially uncorrelated slabs is thus of the form 

TS* 

(f(L,t,T)) = - 1 ^H(t) (47) 

where H (t) is a time dependent Hamaker coefficient whose Laplace transform is given by 

1 f°° / s s \ 

CH(s) = -— uduln(l-A 1 {-i-)A 2 (-i-)exp(-u)) (48) 

07TS Jq \ 11/ 

and the superscript indicates that it is the force for initially uncorrelated slabs. The above expression for the time 
dependent Hamaker coefficient can then be written in terms of the polylogarithmic function Lia(z) = ^vT=i z n /n A to 
give 

CH(s) = JUi 3 (AiH^AaHi)) . (49) 



The static equilibrium value is simply recovered from the pole at s — outside and is given by 

1 

8^r 



H eg - — Li 3 (A!(0)A 2 (0)). (50) 



Eq. (|47p represents an interesting and fundamental result. It states that the non-equilibrium force has the same 
separation dependence as the equilibrium force, but its Hamaker coefficient is time dependent. At least on the non- 
retarded level the dynamic effects do not modify the spatial dependence of the van der Waals force and the force is 
instantaneously long range. This decoupling of the spatial and temporal behavior of the force is in sharp contrast to 
that of critical Casimir force for free scalar fields with Dirichlet or Neumann boundary conditions undergoing model 
A dynamics. In that case the average force exhibits diffusive scaling, behaving as f{L,t) = t~ a g(L/\/i) [a, 0,113, El- 
We finally note that the mapping of the Laplace transform of the time dependent force onto an effective equilibrium 
problem is reminiscent of the results found in [ToL lll| when a similar correspondence occurs for the time dependent 
thermal Casimir force for a free scalar field. 

Let us now investigate some special cases of the above general results. 



A. Short time behavior for Debye type dielectrics 



To begin with we consider the temporal evolution of the Hamaker coefficient in the simplest case where the two 
bounding surfaces are described by identical Debye type dielectric responses and the time evolution is limited to short 
time scales. Assuming that there is only one dipole type in each dielectric slab, we thus have 



8 



By dimensional analysis we see that we can write Ki = —, where Tj is a microscopic polarization relaxation time in 
slab i. This change of notation then yields the familiar Debye formula for the dielectric constant 

£ M= eo + T ^-, (52) 

where e(0) — eo = Ae = \. From this we obtain 

<K) = «" + rr#' < 53 > 

which in turn yields 

Ai U s ) = . (54) 

At short times the temporal behavior of H(t) can be obtained by looking at the large s behavior of CH{s). In this 
limit we have that 

As the Ai (—if) are small we are effectively in the dilute two body limit, this means that the correlations other than 
two body ones set in at later time scales. We thus find 

CH(s) « (ei ' £0 3 )(£2 7 o) (56) 
and thus inverting the Laplace transform we find that at short times 

H{t) « (ei 7" )(£2 - 2 £0) e . (57) 

The above initial growth is quadratic in time, reflecting the need for the polarization fields to become correlated. This 
result can be straightforwardly generalized to several polarization types Vi in each slab and we find that at early times 

H(t) « -^-j Y ^i^. (58) 

From this expression one obviously discerns the pairwise nature of the interaction between the distinct polarization 
types v\ and v 2 in slabs 1 and 2. 

B. Long time behavior for Debye type dielectrics 

Next we consider the long time behavior of the Hamaker coefficient in the case of two identical Debye type materials. 
As we know the initial (and final equilibrium) values of the force we can examine its full temporal evolution by analyzing 
the temporal derivative H(t) = dH/dt. As the initial value of H is zero, standard results on Laplace transforms give 

CH{8) = ^Li 3 (A 1 (-i|)A 2 (-i|)) . (59) 

In this case we have 

Ai(-t£)=A 2 (-»£)= £ ~ e ° 1 with a= ( -^. (60) 

Now using the series representation of Li 3 we can invert the Laplace transform term by term to find 

1 ^ /e-e \ 2n t 2 "- 1 exp(-at) 1 . /(e-e )i 



8?r ^ V e r J n 3 (2n-l)! int PV ' \ e r /' 
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where 

R{u) = f = 4 f ds (ln(«) - ln( S )) C0Sh(5) = 1 * for u _> (62) 

^ n 2 (2n)! Jo s u 2 

n— 1 x ' u 

The asymptotic form for u — > oo was derived by expressing i?(it) in terms of hypergeometric functions. Putting this 
all together yield for large t 

Therefore at late times the asymptotic form of the time dependence of the Hamaker coefficient turns out to be 

^)^ e9 -^(-i2_) 2 exp(-f). (64) 

The apparent divergence in the second term above when e — ¥ eo appears strange at first sight, however we must bare 
in mind that the asymptotic expansion we carried out to obtain this result depended on the variable u = (e— eo)t/reo 
being large. We can numerically verify the validity of the asymptotic expansion by comparison with a direct numerical 
evaluation of Eq. (|6Tj) . It is found to be correct but its realm of validity is for very large values of u of the order of 
50. The asymptotic expansion is thus of limited use and just shows that at very late times the final relaxation to the 
equilibrium Hamaker coefficient is exponential with time scale r* = r/2, interestingly independent of e. 

A particularly non-trivial point about the above calculation is that though the pairwise approximation is valid for 
the equilibrium Hamaker coefficient when Ae = e — eo is small, the pairwise approximation cannot be used to extract 
the temporal behavior of the out off equilibrium Hamaker coefficient. This is because Ae appears multiplied by the 
time t and thus the product of the two eventually must become large. Thus even when the final equilibrium result is 
dominated by pairwise interactions, the dynamical evolution to the equilibrium actually depends crucially on the full 
n-body interactions. 

In the case where the two media are of Debye type but with different dielectric parameters the inversion of the 
Laplace transform of H (s) is more complicated and in general we have not been able to find an analytical expression 
as in the case where both slabs are composed of identical dielectric media. However in the case where the media are 
such that 

£i+eo £2 + £o / cc n 

a = ai = = d2 = , (65) 

neo T- 2 e 

this means that Ai(— i|) and A 2 (— £#) have the same poles (at s = —a% = —02), we obtain 



[ \/(ei ~ £o)(e2 - e )t 



Hit) = JL cM - at)R ^- £o — - j (66) 

Here the late time relaxation to the equilibrium Hamaker coefficient is again exponential with time scale r* = 

l/(a — yj(a — ^-)(a — 7^)), which we see is dependent of the dielectric constant but only through the variable a. 

A sub- variant of this situation, which can be analytically resolved, is the case where one of the systems has a much 
shorter relaxation time than the other, for instance ti <S t\. In the inversion of the Laplace transform we may use 
the approximation 



A, - 



f-d) = 62 ' £ ° « = A 2 (0), (67) 

V 2/ e 2 + e + e T 2 s e 2 + e 

i.e. this is essentially the assumption that the polarization field in the slab 2 instantaneously equilibrates with the 
electric field produced by the polarization field of the slab. In addition the approximation is only valid for times t 
such that t 3> n . Within this approximation we find that 

Hit) = JL exp(-a lt )W ( , (68) 

where 



8irt V eon 



W (u) = > ^ 4-5 ~ for u -> 00, (69) 

v ' ^ n.\n 2 11 2 

n=l 
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in the large u limit. Thus for large t we remain with 

H{t) e ° - Vexp f- 2 t 1+C2 ^ • (70) 

w 47rt 2 VA 2 (ei-e )7 p \ n(e 2 + e )/ V ' 

Therefore even though the relaxation time of the slab 2 is very small, and the overall relaxation time scale is set by 
Ti, i.e. of the relaxation time of slab 1, the relaxation time of the full Hamaker coefficient scales with n but depends 
also on the dielectric properties of slab 2 through the factor r^q^j • 



C. Havriliak-Negami type dielectrics 



It is possible to obtain non-Debye like behavior of the dielectric constant by choosing a suitable distribution of 
polarizability and relaxation times for the associated polarization fields, i.e. by assuming the existence of a distribution 



p(x, T ) = ^2 S ( T - r ») (5 (x - Xv) 



(71) 



such that 



c(cj) = e + / drdx p(x,r) 



X 



1 + iuiT 

A very general phenomenological formula for the dielectric constant is the Havriliak-Negami formula [HI, [24 

Ac 



e(w) = e + 



[1 + (iLUT )<*]P ' 



(72) 



(73) 



with a £ [0, 1] and /3 > 0. The Havriliak-Negami form reduces to the Debye-model in the case where a = p = 1. 
When f3 = 1 it gives the Cole-Cole formula and when a = 1 is gives the Cole-Davidson formula. We should note that 
the Havriliak-Negami dielectric function can be written explicitly as a superposition of individual Debye relaxations 
and thus the study of this functional form within the dynamical formalism presented here is justified. In terms of our 
model this model is composed of polarization fields of the same polarizability \ (= Ae = e — eo) but with different 
relaxation times r. This means that the dielectric function can be written in the form 



e(w) = eo + Ae 



X 



where 



P(r) 



m^sin(/W) 



and 



+ iwr 



= tan 



To y TQ ' 

In the short time limit (corresponding to large s) we find 

2 Ql ' 3l+Q2 ' 32 - 5 AeiAe2 



sin(7ra) 



(^) Q +cos(7ra) 



(74) 



(75) 



CH{s) 



7T6 T Q1 T 02 e Q 



which after inverting the Laplace transform gives for short times 

2ai/3i+a 2 /32-5 £ Q i0i+ Q 2/32 



H(t) 



(76) 



(77) 



where T(z) is the Euler gamma function. We thus see that the exponents a and /3 in the Havriliak-Negami formula 
control the early time growth exponent which depends on the product of the two a/3. We can also verify that this 
general formula agrees with Eq. (1581) in the Debye case where all a and /3 are equal to one. 

The late time decay to the equilibrium Hamaker coefficient between two slabs with Havriliak-Negami dielectric 
functions can be extracted from the small s expansion of CH(s) and we find 



H(t) « H eq - ^Li 2 (A 1 (0)A 2 (0)) ( ^ 

7T V(ei + e )r(-l - ai) V 2t / 



(ea + e )r(-l - aa) V 2t > 



(78) 
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FIG. 1: Color online. Time derivative of the Hamaker constant H(t) for initially uncorrelated Debye slabs made up from 
identical media all with with r = 1 and (i) e = 5eo (solid black line) (ii) e = 5eo (short dashed red line) (iii) e = 20eo (long 
dashed green line) 

The late time decay is thus dominated by the term above with the smaller value of cti- This form of the Hamaker 
coefficient relaxation is very different from the pure Debye exponential relaxation, Eq. (|64l) . and obviously shows a 
long time algebraic tail. We thus see that the time needed to relax to equilibrium in this case can be much longer than 
the microscopic time scales r oi that set the characteristic time of Hamaker coefficient relaxation in the pure Debye 
model. 



VI. NUMERICAL RESULTS 

In this section we will numerically compute the time evolution of the Hamaker coefficient to its equilibrium value 
by numerically inverting the Laplace transform H(s). 

A. Debye dielectrics 

The most complete analytical results we have obtained are for the temporal evolution of the Hamaker constant for 
two Debye dielectrics composed of the same material. We have a closed form expression for the temporal derivative 
H{t) given by Eq. ([tH]). In Fig. ([TJ we have compared the analytical form of Eq. ([tJT|) with the result obtained by 
numerical inversion of the Laplace transform for H(t) for the cases r = 1 and e = 5, 10 and 20. Only the analytical 
curves are shown as the accord with the analytical formula is perfect. The figure shows that the temporal derivative 
increases mots rapidly at short times for the systems of higher dielectric constant as predicted by Eq. (|56[) . However 
the final relaxation is slowest for the systems of lower dielectric constant. 

Now we consider a pure Debye case where each slab is characterized by a single time scale Tj and dielectric constant 
e.;. We examine the case where t\ = 10e and e 2 = 2eo in the following three cases where (i) t\ = 1, r 2 = 1 (ii) T\ = 1, 
r 2 = .1 and (iii) t\ = .1 t 2 = 1 (therefore we are measuring time in the units of the larger of the times T\ and t 2 . 
Note that the equilibrium H eq value of the Hamaker coefficient is the same in all of these cases . In Fig @ we plot 
the functions H(t)/H eq for each cases. We see that in the last two cases, where the shorter time scale r = 0.1 is 
introduced, the approach to the equilibrium value is quicker. However the quickest relaxation occurs when the shorter 
of the two relaxation times is associated with the more dilute dielectric medium, i.e that with the lower dielectric 
constant. 
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FIG. 2: Color online. Time evolution of Hamaker constant normalized by its equilibrium value for two initially uncorrelated 
slabs pure Debye-like dielectric slabs 1 and 2 with e\ — 10eo, £2 = 2eo and relaxation times (i) n = 1, Ti = 1 (solid black line) 
(i) r\ = 1, T2 — 0.1 (short dashed red line) (iii) (i) n = 0.1, Ti = 1 (long dashed green line) 

B. Havriliak-Negami dielectrics 

We now consider the Nevriliak-Negami form for the dielectric constant of two medium. For simplicity we consider 
media of the same type with a fixed Ae = e — eo an d we will use units such that tq = 1. The case where Ae = 3 and 
when (3 = 1 (the Cole-Cole) case is shown in Fig. ([3]) for several values of a. We see that as a decreases from 1 toward 
zero, the relaxation to the final equilibrium Hamaker coefficient (which is the same for all the curves as the static 
dielectric constants are the same) becomes increasingly slow as one would expect from Eq. (|78j) . Also, numerical 
fitting of the late time decay toward 1 in Fig ([3]) is compatible with the analytic prediction of the late time exponents 
given in Eq. (fT5|) . The initial behavior is the inverse, the systems with smallest a have a H(t) which grows faster, in 
accordance with the predictions of Eq. ([77| . An interesting feature of Fig ([3]) is that all the curves cross each other 
at the same, isosbestic-like, point in time, at around t = 0.25. 

In Fig (j4} we show the temporal evolution of H(t) for Cole-Davidson type dielectrics for different values of j3 but 
with the same dielectric constants. As predicted from Eq. (I77|) the systems with the smaller values of (3 show the 
most rapid growth in H{t) at short times but in contrast with the Cole-Cole case these systems also converge most 
rapidly to the equilibrium value, and there is no isosbestic-like point in this case. 

VII. DISCUSSION 

We have formulated a theory describing how the thermal van der Waals force between two dielectric objects evolves 
with time towards its equilibrium value. As an example of a general approach we analyze the dynamics of the thermal 
Casimir or zero Matsubara frequency van der Waals interactions between two slabs brought into close proximity from 
an infinite separation, corresponding to the limit of the non-quantum part of the total van der Waals interaction. 
Despite these simplifications our calculations give a useful intuition about fluctuation induced forces out of equilibrium. 
The particular strength of our approach is that the Laplace transform (with respect to time) of the dynamical force 
can be expressed in terms of equilibrium results available from a wide range of equilibrium methods. 

We show that for the zero frequency Matsubara term in planar geometry the coarse-grained dynamics of the 
interacting dielectric media enters only via the Hamaker coefficient, while the dependence on the spacing between the 
dielectric interfaces remains unchanged and coincides with the equilibrium scaling. This is a fundamental result but 
is limited to non-retarded form of the interaction only. It is however this part of van der Waals interaction that is 
most important in (bio)colloid- and nano-systems. 

The time evolution of the non-retarded van der Waals force between two surfaces of area S at temperature T 
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FIG. 3: Color online. Time evolution of Hamaker consant normalized by its equilibrium value for two identical initially 
uncorrelated Cole-Cole (Nevriliak-Negami with ft — 1) dielectric slabs 1 and 2 with e = 10eo, and characteristic relaxation 
times r = 1 for different values of a: (i) a = 0.8 (solid black line) (ii) a — 0.6 (red short dashed) (iii) a = 0.4 (green long 
dashed) (iv) a = 0.2 (blue dot-dashed). 



separated by L, f(L, t, T), is then given by 

(f(L,t,T)) = -—H(t). (79) 

While the separation into an equilibrium separation scaling and a non-equilibrium Hamaker coefficient appears to be 
universal for non-retarded interactions, the form of the time dependence of the Hamaker coefficient, H(t), is specific 
and pertains to the dielectric response model of the two bounding dielectric surfaces. 

We have shown that for Debye-type dielectric response of the interacting materials with a relaxation time of t, the 
scaling of the Hamaker coefficient for short times is given by 

H(t) ~ C< i 2 , 

while in the asymptotic time regime we obtain 

C If 

t Z T 

The time scale of these non-equilibrium effects in thermal Casimir interactions is thus determined by the (longest) 
dielectric relaxation time of the interacting media. This would render the practical observation of these non-equilibrium 
effects difficult in general, however it may be possible to observe temporal evolution of the force for systems with 
extremely long relaxation times such as polymers and colloids and glassy systems. 

Indeed in the case of the non-Debye-like response that we have studied , i.e. the Havriliak-Negami dielectric response 
function (which is commonly applied to polymeric systems), we obtain a completely different asymptotic behavior of 
the Hamaker coefficient. Instead of an exponential scaling with the relaxation time, a long algebraic tail is obtained 
instead. We derived the scaling form 

in the late time asymptotic regime, where a is one of the scaling exponents in the Havriliak-Negami dielectric response 
function. This long time algebraic tail in the relaxation of the non-equilibrium Hamaker coefficient leads to the 
conclusion that for this particular dielectric model it might be possible to observe long time non-equilibrium effects. 

Variation in the form of the time evolution of the non-equilibrium Hamaker coefficient with the nature of the 
dielectric response of the interacting media makes it possible, or indeed quite probable, that in some experiments 
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FIG. 4: Color online. Time evolution of Hamaker consant normalized by its equilibrium value for two identical initially 
uncorrelated Cole-Davidson (Nevriliak-Negami with a = 1) dielectric slabs 1 and 2 with e = 10eo , and characteristic relaxation 
times r = 1 for different values of a: (i) j3 = 0.8 (solid black line) (ii) f3 = 0.6 (red short dashed) (iii) /3 = 0.4 (green long 
dashed) (iv) (3 = 0.2 (blue dot-dashed). 



where at least indirectly a time-dependent van der Waals interactions are probed, these effects may complicate a 
clear cut interpretation of the experiments. This would be especially true for the tapping mode AFM measurements 
of macromolecular interactions or any other situation involving time varying separation between the interacting 
dielectric interfaces. The AFM tapping mode vibration of the interacting surfaces together with the time evolution of 
the appropriate non-equilibrium Hamaker coefficients would lead to a non-trivial modification of the dependence of the 
magnitude of the van der Waals interaction on the dielectric characteristics of the interacting surfaces. Interestingly 
enough, at least for the non-retarded component of the interaction, the spatial dependence should remain intact. The 
non-equilibrium effects considered in this paper would thus complicate the quantitative interpretation of experiments 
on van der Waals interactions in terms of the dielectric susceptibilities of the interacting interfaces. In order to interpret 
these experiments one would thus possibly need to incorporate the theory presented above to extract quantitative 
measures of time-averaged van der Waals interactions. 

Another possible area of further analysis would be a comparison between these dynamical non-equilibrium van 
der Waals forces with the hydrodynamic drag forces of the Stefan type. These forces have their origin in the flow 
of fluid from between the two interacting surfaces as they are pushed together (25[. The time dependence of the 
non-equilibrium van der Waals interaction, and the specific model studied here, could also be related to, and shed 
light on, the phenomenon of drag induced on moving bodies by Casimir forces [26h29| . 
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